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Abstract. We present a model for the radio outbursts of 
microquasars based on the assumption of quasi-continuous 
jet ejection. The jets are 'lit up' by shock fronts trav- 
eling along the jets during outbursts. The shocks accel- 
erate relativistic particles which emit the observed syn- 
chrotron radiation. The observed comparatively flat decay 
light curves combined with gradually steepening spectral 
slopes are explained by a superposition of the radiation 
of the aging relativistic particle population left behind by 
the shocks. This scenario is the low energy, time-resolved 
equivalent to the internal shock model for GRBs. We show 
that this model predicts energy contents of the radiating 
plasma similar to the plasmon model. At the same time, 
the jet model relaxes the severe requirements on the cen- 
tral source in terms of the rate at which this energy must 
be supplied to the jet. Observations of 'mini-bursts' with 
flat spectral slopes and of infrared emission far from the 
source centre suggest two different states of jet ejections: 
(i) A 'mini-burst' mode with relatively stable jet produc- 
tion and weak radio emission with flat spectra and (ii) 
an outburst mode with strong variations in the jet bulk 
velocities coupled with strong radio emission with steeper 
spectra. We also show that the continuous jets in micro- 
quasars should terminate in strong shocks and possibly 
inflate radio lobes similar to extragalactic jet sources. We 
investigate the possibility of testing the predictions of this 
model with resolved radio observations. Finally, we sug- 
gest that Doppler-shifted X-ray iron lines, and possibly H- 
alpha lines, may be emitted by the jet flow of microquasars 
if thermal instabilities analogous to those in SS433 exist 
in their jets. 

Key words: Stars: individual: GRS 1915+105 - Stars: 
binaries: general - radio continuum: stars - Shock waves 



1. Introduction 

The number of stellar mass black hole candidates known to 
produce jets has increased considerably over recent years. 
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While SS433 was for a long time thought to be a rather 
exotic object, advances in X-ray astronomy have dramat- 
ically increased the number of known galactic X-ray bi- 
naries (van Paradijs 1995). At least nine of these sources 
subsequently showed evidence for the production of rela- 
tivistically moving jets. Mirabel & Rodriguez (1999) give 
an excellent review on the observational and theoretical 
status of these objects, the microquasars. 

The presence of relativistically moving material in 
these sources was discovered during radio observations. 
During times of strongly enhanced radio emission, which 
in the following we will refer to as radio outbursts, the 
emission region can often be resolved into at least two 
components. These components are observed to separate 
in opposite directions over a fews tens of days (i.e. Mirabel 
& Rodriguez 1994). The projected velocity of the compo- 
nent traveling on a trajectory towards the observer can 
exceed the speed of light, if its intrinsic velocity is large 
(e.g. Rybicki & Lightman 1979). In addition to the ap- 
parently superluminal nature of the jet components, their 
propagation has been observed to slow down only in one 
case (XTE J1748-288, Hjellming et al. 1999). This con- 
stant expansion speed led to the interpretation of prac- 
tically ballistic trajectories of discrete plasmon ejections 
as explanation for the observed radio components (e.g. 
Mirabel & Rodriguez 1999). 

Although suggested for the jets of SS433 (Hjellming 
& Johnston 1988), models with quasi-permanent jet pro- 
duction have received little attention in the case of mi- 
croquasars. This is somewhat surprising given the large 
number of similarities they show with jet producing ex- 
tragalactic sources like quasars (hence the name micro- 
quasars) and radio galaxies (e.g. Mirabel & Rodriguez 
1998). In these cases there is little doubt that apart from 
some possible minor intermittency of the jet production 
mechanism for some sources (i.e. Reynolds & Begelman 
1997) the jet flow is practically continuous. 

In this paper we endeavour to close this gap by the 
development of a continuous jet model for microquasars 
to explain the radio outbursts observed in these sources. 
The model is based on the idea of internal shocks in jets 
(Rees 1978) which was successfully applied to Gamma 
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Ray Bursts (GRB) (Rees & Meszaros 1994). In Sect. | we 
briefly review the plasmon model and discuss possible im- 
provements on this within the jet picture. We develop the 
treatment of the relativistic jet flow in Sect. [3] and discuss 
the evolution of the synchrotron emission resulting from 
the internal shock in Sect. [|. The model is then applied 
to the probably best studied radio outburst of any micro- 
quasar, the March 1994 event in GRS 1915+105 (Mirabel 
& Rodriguez 1994) in Sect. ||. The properties of the jet of 
GRS 1915+105 like its energy content are derived from the 
model in Sect. [| In Sect. Qwe consider the implications of 
a continuous jet for the interaction of microquasars with 
their environment. Finally, in Sect. ||we summarise some 
observational consequences of our model which may be 
used to test the model with future observations. 

2. Plasmons or jets? 

Atoyan & Aharonian (1999) developed a model which is 
intended to reconcile the idea of discrete magnetised plas- 
mon ejections during radio outbursts with the observed 
lightcurve and spectral behaviour in the radio waveband. 
They find that a single population of relativistic particles 
accelerated at the time of the ejection of the plasmons 
cannot explain the observations. The energy losses of the 
relativistic particles due to synchrotron radiation would 
lead to a sharp cut-off in the radio spectrum moving to 
lower frequencies as the plasmons expand and travel out- 
wards. This is quite different from the observed rather 
gentle steepening of the spectrum. Atoyan & Aharonian 
(1999) also show that a continuous replenishment of rel- 
ativistic particles to the plasmons alone cannot solve this 
problem because in this case the spectral cut-off moves 
to higher frequencies with time. They therefore postulate 
that the relativistic particles in the plasmons during the 
March 1994 event in GRS 1915+105 were continuously re- 
plenished, presumably by a shock at the side of the plas- 
mons pointing towards the source centre, but also suffered 
energy dependent escape losses. 

To fit the observations the scenario proposed by 
Atoyan & Aharonian (1999) requires that the mean free 
path of the most energetic relativistic particles in the mag- 
netised plasmons is comparable to or exceeds the physical 
dimensions of the plasmons. This implies that these par- 
ticles travel through the plasmons producing synchrotron 
emission but then leave them without scattering once off 
irregularities in the magnetic field or other particles. This 
is difficult to reconcile with the requirement that in order 
to be accelerated to relativistic velocities in the shock re- 
gions the mean free path in these regions must be short 
to ensure many shock crossings. If the accelerating shocks 
are close to the plasmons, this then means that the prop- 
erties of the plasma change dramatically over short dis- 
tances. Moreover, it is not clear in this scenario why the 
synchrotron emission is not completely dominated by the 
contribution of the shocks themselves. 



In this paper we propose a different scenario to explain 
the observed properties of the radio emission of micro- 
quasars during outbursts. This is based on the assump- 
tion that microquasars may produce continuous jets for a 
long time before the actual outburst occurs and may well 
do so permanently (see also Levinson & Blandford 1996a, 
b). The outbursts in our model are then caused by two 
shocks traveling along these continuous jets which accel- 
erate the required relativistic particles in situ. After the 
shock has passed a particular region in one of the jets, this 
region continues to contribute to the total emission until 
the cut-off in the specific spectrum of this region moves 
below the observing frequency. The jet components ob- 
served in microquasars are in general not resolved and so 
the measured flux is the integrated emission from all the 
jet regions passed by the shock which are still emitting 
at the relevant frequency. This implies that the observed 
spectrum is steeper than that of the jet region immediately 
behind the shock where radiation losses are still negligi- 
ble. The variation of the strength of the magnetic field and 
of the number of relativistic particles accelerated by the 
shock along the jet then give rise to a slowly steepening 
radio spectrum. This effect was discussed in the case of 
the radio hot spots of powerful extragalactic radio sources 
by Heavens & Meisenheimer (1987). 

Our investigation is based on the jet model of Bland- 
ford & Rees (1974) and Marscher & Gear (1985). A similar 
approach to explain the synchrotron self-Compton emis- 
sion of extragalactic jets was taken by Ghisellini et al. 
(1985). They develop a numerical scheme to follow the 
evolution of the energy spectrum of the relativistic par- 
ticles downstream of the jet shock taking into account 
radiative as well as adiabatic energy losses. They also in- 
clude synchrotron self-absorption in this calculation. Since 
we are mainly interested in the radio emission of the jets 
of microquasars on rather large scales (~ 10 14 m), we can 
neglect any Compton scattering and absorption effects on 
the energy spectrum. In this case only adiabatic and syn- 
chrotron losses are important and we can use the analytic 
solution for the evolution of the energy spectrum of the 
relativistic electrons derived by Kaiser et al. (1997). 

The underlying physical processes of the model pre- 
sented here are very similar to the internal shock model 
proposed as explanation for GRB (Rees & Meszaros 1994). 
The same scenario has also been invoked to explain the 
X-ray and 7-ray emission of extragalactic jets (Ghisellini 
1999). In the internal shock model the energy of the shock 
traveling along the jet is thought to be supplied by the col- 
lision of fast shells of jet material with slower ones (Rees 
1978). In the case of GRB this energy is released practi- 
cally instantaneously leading to the extremely short du- 
ration of the observed bursts of emission. In extragalactic 
jets the shell collision may take longer but the large dis- 
tance to these objects makes it difficult to separate the 
contributions of multiple collision to the total emission. 
As outlined above, we argue in this paper that the jets 
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of microquasars provide us with the possibility to observe 
the development of internal shocks in jets resolved both 
in space and time. 

3. Dynamics of the jet 

We follow Marscher & Gear (1985) and Ghisellini ct al. 
(1985) in the assumption that all relevant physical quan- 
tities of the jet and the jet material are simple power law 
functions of the unprojected distance from the source cen- 
tre, R. The shock traveling along the jet will compress the 
jet material but we will assume that it does not change the 
behaviour of the physical quantities as a function of the 
distance from the source centre. The radius of the cross 
section of the jet is assumed to follow r oc R ai . For a freely 
expanding, conical jet a\ — 1. The energy density of the 
magnetic field as measured in the rest frame of the jet 
material is then given by u' B ~ u' B (R )(R/ R ) a2 , where 
energy flux conservation would require ci2 = — 2a\. Here 
and in the following dashes denote quantities measured 
in the rest frame of the jet material moving at relativis- 
tic speeds while quantities measured in the frame of the 
observer are undashed. 

The two frames of reference, the rest frame of the ob- 
server and that of the shocked jet material, are defined 
such that the origins of both coincide when the radio out- 
bursts starts, i.e. when the shocks are formed in the centre 
of the source and start traveling outwards. Consider a sec- 
tion of the jet after the shock has passed through it. At 
time t this section is located at 



i?ob = sin9R = sin9v s t/(l ± /3j cos( 



(1) 



in the rest frame of the observer. Here, 9 is the angle of the 
jet to the line of sight and R is the unprojected distance 
of the jet section from the centre of the source, i.e. the 
origin of the observer's rest frame. v s is the deprojected 
velocity of the shock as measured in the rest frame of the 
observer and /3j = Vj/c is the deprojected velocity of the 
shocked jet material in this frame in units of the speed of 
light. The expression in brackets in Eq. ([!]) takes account 
of the Doppler shifted time measurements taken in the 
observer's frame caused by the receding (+) or approach- 
ing (— ) component of the motion of the jet material (e.g. 
Rybicki & Lightman 1979). Transforming t and R to the 
frame comoving with the jet material we find 



t' 
R' 



7j 



V] R 



1 ± /3j cos 9 c 

- \R v jl 

,J V l±/3jCos6» 



(2) 
(3) 



where 7j is the Lorentz factor corresponding to the veloc- 
ity of the shocked jet material, v- y For the origin of the 
comoving rest frame, R' = 0, we recover from Eqs. (|J) 
and (0) the well-known result t' = tS± with the usual 



relativistic Doppler factor S± — (1 ± /3j cos#)] _1 . Since 
the shock moves along the jet, most of the observed emis- 
sion is not produced at the origin of the comoving frame 
but at R' ^ 0. Suppose the section of the jet introduced 
above was passed by the shock at time t' s as measured in 
the frame comoving with the shocked gas. Since this sec- 
tion is at rest in this frame, it is subsequently located at 
R' s = v' s t' s , where v' s is the velocity of the shock as measured 
in the shocked gas' frame. Using this and substituting R 
from Eq. (|) in Eq. §) yields 



t' = tS-l 



(4) 



This expression illustrates the fact that the emission we 
observe at a given time t from different parts of the jet was 
not produced simultaneously in the frame of the shocked 
gas. 

We assume that only those parts of the jet contribute 
to the total emission which have been passed by the shock. 
The emission regions within the jet can therefore be la- 
beled with their 'shock time', tL and are observed at the 
intrinsic time t' given by Eq. (|]). Since t' > t' s > 0, this 
implies t8±c 2 /(c 2 + Vjv' s ) < t' < tS±. For convenience we 
introduce the ratio r{t') = R/R = Jj(v' s t' s + Vjt')/R , 
where we have used Eqs. (||) and (||). 

4. Synchrotron emission of the jet 

4-1- Energy losses of the relativistic electrons 

We assume that the shock passing through the jet material 
accelerates a population of relativistic electrons and/or 
positrons. During the acceleration process and afterwards 
these relativistic particles are subject to energy losses due 
to the approximately adiabatic expansion of the jet ma- 
terial and synchrotron radiation. To determine the exact 
form of the energy spectrum of the relativistic particles 
in a given jet region the kinetic equation including ac- 
celeration and energy terms must be solved. Heavens & 
Meisenheimer (1987) present analytic and numerical solu- 
tions for some simplified cases. They find that the energy 
spectrum follows a power law with a high energy cut-off. 
The cut-off occurs at the energy for which energy gains 
due to shock acceleration balance the synchrotron energy 
losses (e.g. Drury 1983). The cut-off becomes steeper fur- 
ther downstream from the shock. For simplicity we assume 
that the relativistic particles in a given jet region are ini- 
tially accelerated at a time t' s during a short time interval 
dt' s to a power law spectrum with a sharp high energy cut- 
off at 7 max (is). In terms of the number of particles this can 
be expressed by 



N\ ls )d ls = 



K(t' s ) 7 -Pd % dt' s ; 7 s<7max 
; 7s > 7max- 



(5) 
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Here N^(t' s ) is the rate at which relativistic particles are 
accelerated in the jet by the shock at time t' s . The normali- 
sation of this energy spectrum and the position of the high 
energy cut-off depend on the local conditions for diffusion 
in the jet (e.g. Drury 1983). These are not straightfor- 
ward to estimate and we therefore assume for simplicity 
that the initial high energy cut-off of the relativistic par- 
ticles freshly accelerated at time t' s is independent of t' s . 
In the internal shock model for GRB the shock is caused 
by the collision of shells of jet material moving at differ- 
ent velocities. For GRB it is implicitly assumed that the 
collision energy is dissipated very close to instantaneously. 
In the case of microquasars the propagation of the shock 
is resolved in time. For the normalisation of Eq. (^) we 
therefore assume 



(6) 



where R s is the position of the shock at time t' s and 04 a 
model parameter. This implies that the rate at which the 
collisional energy is dissipated and partly conferred to the 
relativistic particles is almost constant for R S /R Q < a<± 
and decreases exponentially at larger distances. The onset 
of the exponential behaviour then signifies the point at 
which almost all of the collisional energy has been dissi- 
pated and the shock starts to weaken significantly. This 
would coincide with the time at which the two colliding 
shells have practically merged into one. Alternatively, in 
only intermittently active sources the exponential decrease 
may be caused by the shock, and therefore the fast shell 
causing the shock, reaching the end of the jet. 

After the passage of the shock the relativistic particles 
continue to loose energy. The rate of change of the Lorentz 
factor of these particles due to the nearly adiabatic expan- 
sion of the jet is given by (e.g. Longair 1981) 



df 



dAV 



3AV dt' 



7 



(7) 



where AV' is the volume of the jet region the particles 
are located in. We assume that the bulk velocity of the 
shocked jet material is constant and this implies that the 
jet is only expanding perpendicular to the jet axis, i.e. 
AV oc (R/R ) 2ai = r 2ai . Changing variables from t' to 
t' then yields 



d-f 2ai 7 
d? ~ 3~ ~ r 

Energy losses due to synchrotron radiation give 



7j«j 



dr'' 3m e c* 7 ^ (jRo) — 1 



(8) 



(9) 



where or is the Thompson cross section and m c the rest 
mass of an electron. By summing Eqs. (||) and @ and 
integrating we find the Lorentz factor 7 at time t' of those 



electrons which had a Lorentz factor 7 S at time t' s (see also 
Kaiser et al. 1997) 



7(f,f. 



with 



//y\-2/3ai 



7sr'(i') 



/(t/)-2/3,x +b 1 (tf,f e h. 



(10) 



hl {t'X) = 4r^ Ub{Ro)Ro [ T '(t'r 3 - r'(Q a °] (11) 



3a 3 m c 7j«j 



and 



&3 = 1 - 2 ( a 2 



(12) 



The number of relativistic particles with a Lorentz factor 
in the range 7 to 7 + dj in the jet region overtaken by the 
jet shock at time t' s is therefore given by 

K(Kh- p b 2 

N'h)dl = { XT l {t')- 2/5ai d 1 dt' s ;7<7ma*(0 (13) 

;7>7max(i')> 



with 



p-2 



r'(t')- 2/3ai -htfXft r'^yl^v-V. (14) 



Note here that the high energy cut-off, 7 ma x(£')j a l so 
evolves according to Eq. (|l0|). 

4-2. Synchrotron emission 

The synchrotron emission of the relativistic particles at 
time if in the region of the jet overtaken by the shock at 
time t' s is given by 



dP!,, = 



-a T cu B {t', t' s h 2 <P(v', j)N'{j)dj, 



(15) 



where $(^',7) is the synchrotron emission spectrum of 
a single electron with Lorentz factor 7. Here we assume 
that the magnetic field in the jet is tangled on scales 
smaller than the radius of the jet. This then implies that 
$(^', 7) = F(j)/v c , where F is one of the synchrotron in- 
tegrals normalised to give J x F(-f)dj = 1 (e.g. Shu 1991). 
To get the total emission of the jet behind the jet shock 
we have to sum the contributions of all the regions labeled 
with their shock times t' s within the jet. From Eq. ( p^ ) we 
see that this implies integrating Eq. (JlJ) over t' s . This in- 
tegration must be performed numerically. Finally, to com- 
pare the model results with the observations we have to 
transform to the rest frame of the observer, P v = Pl,8±. 
Note here that the luminosity inferred from the observa- 
tions at frequency v was emitted in the gas rest frame at 
a frequency v' = v/8±. 

As mentioned above, the steepening of the radio spec- 
trum of the superluminal jet components observed dur- 
ing outbursts of microquasars is caused by the integrated 
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emission from an extended region of the jet. The further 
away from the shock the emission is created in the jet, the 
lower the cut-off in the energy spectrum of the relativistic 
electrons will be. However, if the jet region contributing 
to the total emission is not too large the properties of the 
energy spectrum will not change dramatically within this 
region. In this case we may estimate the approximate loca- 
tion of the break in the radio spectrum beyond which the 
spectrum steepens significantly. Most of the energy lost by 
relativistic particles is radiated at their critical frequency, 
v c = (3/2) ^l7 2 , where z/l is the Larmor frequency. Defin- 
ing the break frequency, z^b, as the critical frequency of 
the most energetic particles just behind the jet shock we 
get 



3qy/2fi u' B 

47TTO C 



7max i 



(16) 



where q is the elementary charge and fj, Q is the magnetic 
permeability of the vacuum. The steepening of the ra- 
dio spectrum of the observed outbursts of microquasars 
strongly suggest that the observing frequency, v, is close 
to the break frequency, i.e. v ~ v[J)±. Therefore, if most 
of the observed emission comes from the region just be- 
hind the shock, we expect from Eq. ( |l6| ) that 7max oc 
Ug(i? ) -1 / 4 . This implies that j max and u' B (R ) are not 
independent parameters of the model but that they are 
correlated. 




Log(t/days) 

Fig. 1. Flux density measurements and model fits of the 
March 19 radio outburst of GRS 1915+105 at 8.4 GHz. 
Crosses and diamonds: Approaching jet component dur- 
ing first and second VLA campaign respectively, stars 
and triangles: receding jet component, same observations. 
Squares: Total radio flux at 8.4 GHz multiplied by a factor 

5 for clarity. Solid lines: Predictions of the fiducial model 
fitted excluding flux measurements at 3.3 GHz for (from 
top to bottom) the total flux, the approaching and the 
receding jet component. Dashed lines: like solid lines but 
using all data for the model fit. Observations from Mirabel 

6 Rodriguez (1994) and Rodriguez & Mirabel (1999). 



5. Application to GRS 1915+105 

The number of parameters in the model outlined above is 
large. In order to reduce this number we assume that the 
jets in microquasars are freely expanding, i.e. a± = 1, and 
that the flux of magnetic energy through the jet is con- 
served, i.e. a 2 = —2. Since we assume the bulk velocity of 
the jet material to be constant, this implies that the ratio 
of the kinetic energy and the energy of the magnetic field 
is constant as well. Furthermore, we impose symmetry be- 
tween the approaching and receding sides of the source in 
the sense that the model parameters describing the jet are 
the same on both sides. This may be a poor assumption as 
the jets of GRO J1655-40 are observed to be asymmetric 
(Hjellming & Rupen 1995). The model then depends on 
five free parameters: The e-folding distance of the number 
of relativistic particles within the jet accelerated by the 
shock, 0,4, the bulk velocity of the shocked jet material, uj, 
the maximum Lorentz factor up to which relativistic par- 
ticles are initially accelerated, 7 maxj the slope of the initial 
power law energy spectrum of these particles, p, and the 
energy density of the magnetic field at R , u' B (R ). The 
acceleration rate of relativistic particles at the normalisa- 
tion radius, N^(R ), is in principle also a free parameter. 
However, from Eq. [l5|) we note that it is only a multiplica- 
tive factor in the calculation of the total radio emission of 
the jet. We therefore use it to normalise the model in such 
a way that for a given set of model parameters N^(R a ) 



is such that the difference between the model predictions 
and the observational data is smallest for this given set of 
parameters. 

Many radio outbursts of a number of microquasars 
have been observed. But to constrain the model param- 
eters in a meaningful way, we would ideally need radio 
observations at two or more frequencies which clearly re- 
solve the approaching and the receding jet component. 
Furthermore, the resolution should be sufficient to decide 
whether one of these jet components consists of multi- 
ple subcomponents, i.e. multiple shocks, the emission of 
which may be blended in observations of lower resolution. 
To date there are very few simultaneous multi-frequency 
observations which come even close to this ideal situation. 
The best studied radio outburst of any microquasar is 
still that of March 19th 1994 of GRS 1915+105 (Mirabel 
& Rodriguez 1994). This is also the outburst studied by 
Atoyan & Aharonian (1999). To test our model we will use 
the comparatively large data base accumulated during this 
event. 

5.1. The observations 

The radio outburst of GRS 1915+105 which occurred in 
March 1994 was one of the strongest recorded for this ob- 
ject. The flux density at 1.4 GHz exceeded 1 Jy which is 
at least ten times higher than the radio flux in quiescence 
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(Rodriguez et al. 1995). The outburst was observed with 
the VLA in A-array at 8.4 GHz during 7 epochs covering 
almost 42 days. Except for the first of these, the two jet 
components were resolved at this frequency (Mirabel & 
Rodriguez 1994). For the first unresolved and one further 
epoch measurements with the VLA are also available at 
4.9 GHz and 15 GHz. In addition, GRS 1915+105 was 
monitored during this time by the Nancay telescope at 
1.4 GHz and 3.3 GHz (Rodriguez et al. 1995). There are 
24 flux measurements at each frequency but the source 
is unresolved at these frequencies. After a reconfigura- 
tion of the VLA four more observations of GRS 1915+105 
of lower resolution were obtained at 8.4 GHz in B-array 
(Rodriguez & Mirabel 1999). These measurements cover 
roughly another 28 days but there is a gap of about 40 
days between the end of the A-array observations and the 
start of the B-array campaign. 

From the resolved VLA observations at 8.4 GHz 
Mirabel & Rodriguez (1994) determined a velocity for the 
jet components of 0.92 c and an angle of the jets to the line 
of sight of 70° . This assumes that the approaching and the 
receding component travel at the same velocity in opposite 
directions. Furthermore, it is assumed that GRS 1915+105 
is located 12.5 kpc away from us. Fender et al. (1999) ob- 
served another radio outburst of GRS 1915+105 in Oc- 
tober 1997 and found a higher intrinsic velocity which is 
inconsistent with a distance of 12.5 kpc. They argue that 
the most likely distance for this object is 11 kpc which 
then implies that the velocity of the jet components in 
March 1994 was 0.86 c and the angle of the jets to the 
line of sight is 68°. In the following we will adopt these 
later values. Extrapolating back the trajectories of the two 
jet components Mirabel & Rodriguez (1994) find that the 
outburst started at 20 hours on March 19. 

Figs. and || show all available flux density measure- 
ments as a function of time. Rodriguez & Mirabel (1999) 
note that another outburst of GRS 1915+105 occurred on 
April 21. The jet components of this new outburst are 
clearly visible as an unresolved emission peak coincident 
with the source centre in the VLA radio map of epoch 6. 
The approaching component of this new outburst is also 
distinctly visible on the map of the following observing 
epoch while the receding component is probably blended 
with that of the previous outburst of March 19. The sig- 
nature of this later outburst as a sudden increase of the 
radio flux is not very distinct at 1.4 GHz and 3.3 GHz. 
Even at 8.4 GHz the situation in terms of the total flux 
density is somewhat unclear. However, the outburst can 
be easily identified as a separate event from the one of 
March 19 because the VLA maps reveal an emission peak 
distinct from those of the earlier outburst. The later ob- 
servation epochs at the VLA with lower resolution detect 
the jet components of the April 21 burst while the compo- 
nents of the March 19 event were not detected (Rodriguez 
& Mirabel 1999). This is rather puzzling as the extrap- 
olation of their lightcurves from the earlier observations 
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Fig. 2. Flux density measurements of the March 19 ra- 
dio outburst of GRS 1915+105 at 1.4 GHz and 3.3 GHz. 
Crosses: Total radio flux at 1.4 GHz, stars: Total radio flux 
at 3.3 GHz divided by a factor 5 for clarity. Model predic- 
tions as in Fig. [I]. Observational data from Rodriguez et 
al. (1995). 

indicate that they should still have been visible during 
these later observations. 

5.1.1. The blending of outbursts 

A close examination of the total radio flux at 8.4 GHz re- 
veals another sudden increase around observation epoch 
4 (April 9, Fig. 0). This increase is also detected at 1.4 
GHz and 3.3 GHz by the Nancay telescope on April 5 
(Rodriguez et al. 1995, Fig. ||). Atoyan & Aharonian 
(1999) point out that this may be caused by a sudden 
additional injection of fresh relativistic electrons in the 
'blobs' of gas they consider in their model. An alternative 
interpretation, which we will adopt here is that another, 
smaller outburst occurred shortly after April 4. In this 
case, the radio emission caused by the propagation of a 
new shock down the jet is most likely blended with that 
of the earlier event of March 19 because of the limited 
resolution of the observations. If the smaller outburst oc- 
curred at April 5, 0:0 hours, and assuming that the appar- 
ent shock velocity of this outburst is equal to that of the 
March 19 event, i.e. 17.5 mas day -1 (Rodriguez & Mirabel 
1999), then the distance of the shock from the source cen- 
tre on the approaching jet side at observing epoch 4 would 
be roughly 0.1". The position of the emission peak on 
April 9 is given as 0.36" from the source centre (Mirabel 
& Rodriguez 1994). The resolution of the VLA in A-array 
at 8.4 GHz is ~0.3" and this means that a secondary peak 
caused by a later, somewhat weaker outburst as proposed 
here could not be detected as an individual structure. The 
situation on the receding jet side with its lower expansion 
velocities is even worse. However, the radio emission of 
the additional outburst would contribute to the total ra- 
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dio flux of the source and we believe that this has been 
detected here. 

In summary, we will assume that the observations out- 
lined above cover three separate radio outbursts of GRS 
1915+105. The first and strongest occurred on March 19. 
The jet components of this burst are detected until the 
end of the first set of VLA observations (April 30). By 
the time of the second VLA campaign (June 13) they 
had vanished although the extrapolation of their earlier 
lightcurves suggested that they should still be observable. 
The second much weaker outburst occurred shortly after 
April 4, since the Nancay data show a sudden increase in 
radio flux on April 5 but not on April 4. The radio emission 
of this event is most likely blended with that of the first 
outburst and there is no sign of this second burst in the 
second set of VLA observations. The third outburst finally 
started on April 21 and was intermediate in strength. The 
radio emission caused by this event is clearly detected in 
the source centre on April 23 and the approaching jet com- 
ponent can be seen in the VLA map obtained on April 30. 
Both jet components are clearly detected in all four VLA 
observing epochs in June and July. 



5.1.2. Data used in the modeling 

In the continuous jet model for microquasars outlined 
above, only one shock is thought to travel outwards in 
each jet. Because of this, only observational data from ob- 
serving epochs during which we can be sure that there was 
only one shock per jet contributing to the radio emission 
can be used in constraining the model parameters. From 
the above discussion it is clear that we can only use the 
three unresolved VLA measurements at 4.9 GHz, 8.4 GHz 
and 14.9 GHz from March 24 and the two following re- 
solved observations at 8.4 GHz. These later observation 
provide us with separate flux measurements for the ap- 
proaching and receding jet components. We also use the 
measurement of the receding jet component of April 9, 
since it seems likely that all the flux of the second outburst 
was attributed to the approaching component by Mirabel 
& Rodriguez (1994). Alternatively, this measurement can 
be taken as an upper limit. Of the Nancay observations 
we use all flux measurements starting March 24 through 
to April 4. There are eight observations during this time 
at 1.4 GHz and 3.3 GHz. For all 24 measurements used 
to constrain the model we assume the conservative error 
of 46% suggested by Rodriguez & Mirabel (1999) as op- 
posed to the original error of 5% quoted by Rodriguez et 
al. (1995). In practice we did not use the eight 3.3 GHz 
data points as their inclusion led to significantly worse 
fits of the model to the observational data. See the next 
section for a discussion of this point. 



5.2. Constraining the model parameters 

We use the model outlined above to calculate the expected 
radio flux at the times GRS 1915+105 was observed dur- 
ing the outburst starting March 19. The 'goodness of fit' 
of the model for a given set of free parameters to the ob- 
servations was assessed by calculating the sum of the \ 2 ~ 
differences at the times the source was observed. The best- 
fitting model was then found by minimising this x 2 -value 
using a 4-dimcnsional downhill simplex method (Press et 
al. 1992). The remaining fifth model parameter, the en- 
ergy density of the magnetic field at the normalisation 
radius, u' B (R ), was set 'by hand' to five different values. 
The normalisation radius, R , was set to 1.1 • 10 14 m which 
corresponds to the unprojected distance of the jet shock 
from the source centre at the time of the first VLA obser- 
vation, i.e. March 24. 

The results of the model fits are summarised in Ta- 
ble [l]. The model fits the observational data equally well 
for all five adopted values for the strength of the mag- 
netic field. As expected from Eq. (|l6|) we find that the 
maximum Lorentz factor up to which relativistic particles 
are accelerated correlates strongly with the value of the 
energy density of the magnetic field. The values found for 
7max in the model fits follow almost exactly a u' B (R o )~ 0/25 
law. This shows that the model presented here cannot be 
used to constrain the strength of the magnetic field within 
the jet independently of the maximum energy of the rel- 
ativistic particles. To proceed we adopt in the following 
u' B (R ) = 4.2 ■ 10 ~ 5 as our fiducial model. This corre- 
sponds to the equipartition value of the magnetic field of 
10~ 5 T (0.1 G) found by Atoyan & Aharonian (1999) for 
the chosen value of R . Note however, that the energy 
density of the magnetic field and that of the relativistic 
particles does not stay in equipartition for all times in our 
model. 

Figs. [j] and || show the model predictions of the fiducial 
model compared to the observational data. Also shown are 
the predictions of the model for the best-fitting parame- 
ters when using also the 3.3 GHz data. If we use all avail- 
able data, the model predictions at the higher observing 
frequencies are rather low at early times. For the VLA 
measurements at 4.9 GHz and 14.9 GHz (not shown in 
Figs. [I] and ||) on March 24 the model predicts flux den- 
sities of 759 mJy and 385 mJy respectively. This is much 
lower than the measured flux densities of 887 mJy and 514 
mJy at these frequencies. A closer inspection of Fig. || also 
shows that the measurement at 4.9 GHz actually exceeds 
those taken at 3.3 GHz at comparable times, which is un- 
likely to be real. Table [l] also shows that the fit obtained 
including the 3.3 GHz data is much worse in terms of the 
reduced x 2 -vahies than that excluding them. Also the flux 
densities predicted by our fiducial model which excludes 
the 3.3 GHz data, 926 mJy at 4.9 Ghz and 499 mJy at 
14.9 Ghz, are much closer to the observations. Finally, 
the predictions of this model at 3.3 GHz fit the observa- 
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Table 1. Results of the model fits. The first five fits are obtained without the 3.3 GHz data points, iig is the energy 
density of the magnetic field, 04 is the radius beyond which the acceleration rate of relativistic particles decreases 
exponentially in units of R (see Eq. ||), wj is the bulk velocity of the shocked jet material, 7 ma x is the initial high 
energy cut-off of the energy spectrum of the relativistic particles, p is the initial slope of this spectrum and is its 
normalisation. 
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tions well at this frequency apart from the early observing 
epochs (see Fig. ||). We therefore believe that the model 
parameters found using our fiducial model and excluding 
the 3.3 GHz data are more reliable than those found when 
including these additional measurements. 

To estimate the expected error of the model parame- 
ters of our fiducial model we calculated the x 2 -value for a 
large set of combinations of the 4 free model parameters. 
The uncertainties quoted in Table ^ are 1-cr errors cor- 
responding to those parameter ranges for which \ 2 < 1- 
Note that the uncertainties of the model parameters, par- 
ticularly those of 04 and p, are large while the light curves 
predicted by the model pass the data points well within 
the error bars of the flux measurements (see Figs, [l] and 
|2|). This suggests that the quoted errors of the observed 
fluxes, at least for the VLA data points, are too conserva- 
tive which also results in an overestimation of the uncer- 
tainties of the model parameters. 

5.3. Comparison with the later observations 

The exponential function which describes the change in 
the acceleration rate of relativistic particles as a function 
of the position of the shock in the jet, Eq. (^J), implies 
that the radio flux caused by the first outburst on March 
19 decreases quickly once R s /R exceeds 04. This effect 
is clearly visible in Fig. [I] at 8.4 GHz. The model pre- 
dicts that without the additional blended radio emission 
caused by the second and third outbursts around April 5 
and April 21 the two jet components would have faded 
much more rapidly than is observed. This effect is less 
pronounced at lower frequencies (see Fig. ||). However, 
even for these the predicted and the observed light curves 
steepen somewhat roughly 15 days after the start of the 
first outburst. The continued steepening of the lightcurves 
of the two jet components at 8.4 GHz can also explain why 
they were not detected during the second observing cam- 
paign at the VLA (see Fig. [l]) . 



We tried replacing the exponential in Eq. (|6j) with a 
simple power law. The data can be adequately fitted with 
this modified model as well. However, this change in the 
temporal behaviour of N^(R Q ) also leads to a much in- 
creased flux at low observing frequencies at later times. 
Using this modified model we found a flux at 1.4 GHz 40 
days after the start of the first outburst exceeding the ob- 
servations by a factor of at least 1.3 even without consid- 
ering the possible contribution from later outbursts. This 
supports the picture of two colliding shells of jet mate- 
rial of finite width causing the internal shock. The rate at 
which energy is dissipated is roughly constant during the 
collision and decreases rapidly once the two shells have 
merged. 

Comparing the model lightcurves with the observa- 
tional data the signature of the second outburst starting 
around April 5, about 16 days after the start of the first 
outburst, can clearly be detected at 8.4 GHz. The increase 
in the radio emission caused by the second event is less 
dramatic at 3.3 GHz and 1.4 GHz but can still be seen 
in Fig. §. The third outburst of April 21, 32 days after 
the first burst, is seen as excess emission at 8.4 GHz and 
3.3 GHz but is less obvious at 1.4 GHz. We note that in 
general the smooth lightcurve predicted by our model fits 
the VLA observations much better than the flux measure- 
ments at 1.4 GHz and 3.3 GHz taken with the Nancay 
telescope. This may imply larger errors for the low fre- 
quency data which hide to some extent the signatures of 
the second and third outburst which are much weaker than 
the first. 

6. Properties of GRS 1915+105 

6.1. Energetics 

Using the model parameters of our fiducial model we now 
derive some of the physical properties of the jets of GRS 
1915+105 during the outburst of March 19. The rate at 
which energy is transfered by the jet shock to the rela- 
tivistic particle population at time t' s is 
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E> el = N^R )e- R °" R °^m c c 2 / 7s ^( 7s - 1) d % 
= io 29 e -^/(^)w, (17) 

where we have assumed that only electrons and/or 
positrons are accelerated and that the initial energy spec- 
trum of the relativistic particles extends down to 7 m i n = 1. 

The rate at which energy is transported in the form of 
magnetic fields can be estimated by 



E' B 



i ^) <+< km 



1.8-10 



2N 



degrees 



W, 



(18) 



where 8 is the opening angle of the conical jet and v' s is the 
speed of the shock in the frame of the shocked jet material. 
For our fiducial model v' s = 0.53 c. Fender et al. (1999) find 
that during another outburst of GRS 1915+105 in 1997 9 
was smaller than 8° . This would then imply an upper limit 
to E' B of 1.2 ■ 10 30 W. Note that the strength of the mag- 
netic field in the jet after the passage of the shock is used 
here. This estimate does not imply that the unshocked jet 
material carries a magnetic field of this strength. Some or 
all of the magnetic field may be generated in the shock 
itself. 

Finally, we can derive an lower limit for the bulk ki- 
netic energy transported by the jet material. We know the 
number of relativistic light particles in the jet and so 



E kin > (7j - l)m c c 2 N' {R )e- R sK R ° a ^ / 7-^7. 

Jfmio 

= 6.3 • io 26 e ^^/(«oa4) w ( 19 ) 

This is only a strict lower limit, since we do not know 
whether the jets also contain thermal material and/or pro- 
tons. In the case that there is one proton for each relativis- 
tic electron we find that the numerical constant in Eq. ( |l9| ) 
increases to 1.2- 10 30 . Note that this then is identical to the 
energy carried in the form of magnetic fields for 9^8°. 

The estimates for the energy transported along the 
jet in various forms presented above are lower by about a 
factor 10 than the estimates of Fender et al. (1999) for the 
weaker outburst in 1997. However, it should be noted that 
their estimates are based on the assumption that the radio 
emission is caused by two 'blobs' of relativistic plasma 
which were ejected by the central source within about 12 
hours. The continuous jet model presented here requires 
that the estimated energy supply to the jet is sustained 
by the central source for at least 42 days; the length of 
the first observing campaign. This means that the total 
amount of energy produced by GRS 1915+105 is predicted 



by our model to be at least an order of magnitude greater 
during the March 1994 outburst than it was in the case of 
discrete ejections assumed for the September 1997 event. 

These estimates illustrate that a continuous jet model 
cannot decrease the total amount of energy needed for a 
given radio outburst but the rate at which this energy is 
produced is much lower than in a model assuming dis- 
crete ejection events. This is the case because much, if not 
most, of the energy needed to produce the radio emission 
observed is 'stored' in the material of the continuous jet. 
This material was ejected by the central source during 
comparatively long period well before the process which 
led to the formation of the jet shock took place. Only the 
acceleration of relativistic particles at the jet shock then 
'lights up' the jet and we are able to detect it. 

6.2. Self- absorption 

Since all jet properties are assumed to scale with distance 
from the source centre in our conical jet, it is clear that 
at some early time in the outburst the jet material was 
opaque for radio emission because of synchrotron self- 
absorption. The absorption coefficient in the rest frame 
of the emitting gas is given by (e.g. Longair 1981) 



18,\P 



X' u , = 3.354 • 10~ 9 (3.54 • 10 



where in our notation 



4K(Ro) (m e c 
nv' s (R s 9) 2 



2\P- 



- P -4)/2 m -l 



-R s /(R oai ) 



(20) 



(21) 



and b(p) is of order unity. We only consider the region just 
behind the jet shock where the energy distribution of the 
relativistic particles is completely described by a power 
law of exponent p. For a photon emitted at the centre of 
the jet the optical depth in the radial direction is then 
r = R s 9x' L ,</2- For our fiducial model we then find that 
the jet material becomes transparent at 8.4 GHz roughly 
2 hours after the start of the outburst when the shock has 
reached a distance of ~ 2 • 10 12 m from the source centre. 

Mirabel et al. (1998) find that for the much weaker 
'mini-bursts' of GRS 1915+105 the jets become transpar- 
ent about 30 minutes after the start of the burst. Bearing 
in mind that the mini-bursts may be quite different in 
their properties compared to the major outburst consid- 
ered here, our value is therefore in good agreement with 
their findings. 

6.3. Infrared emission 

Several groups have reported the detection of infrared 
emission from GRS 1915+105 (i.e. Sams et al. 1996, 
Mirabel et al. 1996, 1998). In the case of the mini-bursts 
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simultaneous flux measurements at radio frequencies and 
in the K-band are available at times of about 10 to 20 
minutes after the start of the bursts (Mirabel et al. 1996). 
Because of the uncertainties in the dust corrections in the 
K-band towards GRS 1915+105 it is difficult to estimate 
the spectral behaviour from radio to infrared wavelengths. 
However, for the mini-bursts flat spectra, a — 0, regard- 
less of the exact magnitude of extinction are observed very 
early during the bursts (Mirabel et al. 1996). The slope 
of the initial energy distribution of the relativistic parti- 
cles in combination with the rather low high-energy cut-off 
of this distribution we found for our fiducial model is in- 
consistent with such flat emission spectra. However, this 
model does predict an unobscured, optically thin infrared 
flux of about 14 mJy in the K-band for a time about 15 
minutes after the start of the burst. This may be enough 
to be detected in future observations of large outbursts. 
The timing requirements for such an observations are how- 
ever difficult to meet, since the predicted infrared flux very 
quickly becomes undetectable at only slightly later times. 

More puzzling is the detection of a resolved jet compo- 
nent with K-band flux of at least 1.8 mJy about 0.3" away 
from the centre of GRS 1915+105 by Sams et al. (1996). 
The shock on the approaching side of our jet model would 
need 24 days to reach such a large distance from the source 
centre. By this time our fiducial model predicts no syn- 
chrotron emission in the K-band at all. We have estimated 
whether this infrared emission may be caused by radio 
photons which are inverse Compton scattered to such high 
frequencies within the jet plasma. However, we find that 
this cannot explain the observations since the density of 
the relativistic particles in the jet in our model is orders 
of magnitude too low. 

The observation of K-band emission far away from the 
core of GRS 1915+105 and the flat spectral indices of the 
mini-bursts suggest that two different types of outbursts 
may occur in the jets of this source. The strong radio 
bursts like the one of March 1994 are caused by jet shocks 
which produce large numbers of relativistic particles with 
a steep energy distribution. The weaker mini-bursts in- 
volve shocks which accelerate less particles but produce 
a flatter energy distribution which may also extend to 
higher energies than in the stronger bursts. The 'mini- 
burst mode' may correspond to a phase of relative stable 
jet production with only small variations in the bulk ve- 
locity of the jet material. Such flat spectra extending to 
millimeter wavelengths, possibly coupled with the contin- 
uous ejection of a jet, have been observed in Cygnus X-l 
(Fender et al. 2000). The strong radio outbursts then prob- 
ably mark phases of more violent changes in the central 
jet production mechanisms. Fender (1999) points out that 
this proposed behaviour may also be reflected in the X-ray 
signature of the accretion disk. In any case, other sources 
of infrared emission in the close vicinity of the jets like 
dust illuminated by the disk and/or the jet may further 
complicate the situation (Mirabel et al. 1996). To test the 



validity of the proposed scenario resolved observations of 
outbursts of GRS 1915+105 and other galactic jet sources 
from radio to infrared frequencies would be necessary. 

7. The end of the jet 

The energy transported by the jets of microquasars is 
enormous. This energy will be continuously deposited at 
the end of the jets and may lead to significant radiation 
from this region depending on how it is dissipated. In the 
following we investigate the fate of the energy transported 
by the jets of microquasars as predicted by our model. The 
discussion is based on the work by Leahy (1991). 

7.1. Momentum balance 

In order for the jets to expand they have to accelerate the 
surrounding ISM and push it aside. The velocity of the 
contact surface between the front end of the jet and the 
ISM, v c , is given by balancing the momentum or 'thrust' 
of the jet material with the ram pressure of the receding 
ISM 



= 1 



T,Mf 



\ 




T.Mf 



(22,) 



where Tj and T c are the adiabatic indices of the jet ma- 
terial and the ISM respectively, Mj is the internal Mach 
number of the jet flow, M c is the Mach number of the con- 
tact surface with respect to the sound speed in the ISM 
and rj = pj p c . Here p- } is the mass density of the jet mate- 
rial while p c is the density of the ISM. This expression is 
strictly valid only for non-relativistic jet velocities. How- 
ever, since the bulk velocity of the shocked jet material, wj, 
is only mildly relativistic in our fiducial model, 7j = 1.3, 
we take Eq. (|2|) to be a good approximation. Note that 
the velocity of the jet material in front of the shock is even 
lower than v\. 



For v c 



the jet material does not decelerate 



strongly at the end of the jet. This implies that little of the 
kinetic energy transported by the jet is dissipated. Even 
for large internal Mach numbers it is then unlikely that 
a strong shock will develop in the jet flow close to the 
contact surface. This occurs when the jet is overdense, i.e. 
rj ^> 1 and so the jet flow is close to being ballistic. For 
underdense jets, r\ < 1, the ratio v c /vj can become consid- 
erably smaller than 1. In this case a strong deceleration of 
the jet ensues and much of its kinetic energy is dissipated. 
For Mj ^> 1 a strong shock will form and can act as a site 
of efficient acceleration of relativistic particles. Examples 
for this are the powerful extragalactic radio sources of type 
FRII (Fanaroff & Riley 1974) with their very bright radio 
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hot spots at the end of their jets. The diffuse radio lobes 
enveloping their jets are the remains of the shocked jet 
material left behind by the advancing contact surface. In 
the transonic regime, A/j ~ 1, only weak shocks may form 
at the jet end and particle acceleration is less efficient. The 
less powerful jets of FRI objects fall in this class. 



7.2. Application to our fiducial model 

In the model developed in the previous sections we have 
assumed the jets of microquasars to be conical with a con- 
stant opening angle. This implies n = r] (R/ R Q )~ 2 and, 
because of the adiabatic expansion of the jet material, 
Mj = M i (R )(R/R ) r '- 1 . The bulk velocity of the jet 
material is high in our fiducial model and unless the jet 
material is very hot (Tj(i? ) > 10 12 K in the case of a 
proton-electron jet) the internal Mach number of the jet 
flow will always greatly exceed 1. Since r\ is a strongly 
decreasing function of R, we expect from Eq. ( p2| ) that 
the ratio v c /vj will always fall significantly below unity 
for large values of R. This means that the jets of micro- 
quasars should end eventually in strong shocks which may 
be detectable in the radio. In the source XTE .11748-288 a 
region of bright radio emission was observed to slow down 
and brighten at the same time some distance from the 
centre of the source (Hjellming et al. 1999). In our model 
this is interpreted as an internal shock reaching the end 
of the jet where the termination shock further boosts the 
relativistic particle population which was pre-accelerated 
by the internal shock. After passing through the termina- 
tion shock the jet material may inflate a radio lobe very 
similar to extragalactic FRII objects if v c /vj <C 1 (see 
also Levinson & Blandford 1996a, b). It has been sug- 
gested that the diffuse radio emission region W50 around 
SS433 is the radio lobe inflated by the jets of this source 
(Begelman et al. 1980). Other radio lobes were detected 
around IE 1740.7-2942 (Mirabel et al. 1992), GRS 1758- 
258 (Rodriguez et al. 1992) and possibly GRO J1655-40 
(Hunstead et al. 1997), but not in the vicinity of GRS 
1915+105 (Rodriguez & Mirabel 1998). The absence of 
a radio lobe in GRS 1915+105 may indicate that the jet 
in this source is relatively young and has not yet reached 
the point at which it becomes underdense with respect to 
the ISM. In the following we estimate the distance out 
to which the jets in this source may travel without the 
formation of a strong termination shock. 

A lower limit for 7/ can be derived from our fiducial 
model assuming that the jets consist only of the relativistic 
particles responsible for the synchrotron emission plus the 
particles needed for charge neutrality. Thus 
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Fig. 3. The advance speed of the termination shock at the 
end of a hypersonic, conical jet in units of the bulk velocity 
of the jet material. Solid line: The case of maximum mass 
transport rate. Dashed line: Minimum mass transport rate 
of a proton-electron jet. Dot-dashed line: Minimum mass 
transport rate for a pair plasma jet. 

the rate at which mass is ejected along the jet can not ex- 
ceed the mass accretion rate within the disk powering the 
jet. Fits to the X-ray spectrum of GRS 1915+105 suggest 
an accretion rate of order m ~ 10 15 kg s _1 (Belloni et al. 
1997). We then find 
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where mj is the mass of the average particle in the jet. An 
upper limit for 77 can be derived from the assumption that 



Note that this upper limit does not depend on the nature 
of the jet material. Using Eqs. ( p3| ) and ( p4| ) and assuming 
Pc ~ TOplO 6 kg m~ 3 , corresponding to a particle density 
of 1 cm~ 3 , we find 56 < r/ < 1300 for a proton-electron 
jet and 0.03 < i] Q < 1300 for a pair plasma jet. The lower 
limit for the pair plasma jet assumes that the pairs are 
cold. Because of pair annihilation it is unlikely that the 
material of a pair plasma jet is cold (e.g. Gliozzi et al. 
1999) and so this lower limit is used here for illustrative 
purposes only. Relativistic thermal motion of the pairs 
would raise this lower limit. 

For the reasonable assumptions M c 3> 1 and Mj 3> 1 
Fig. H shows the ratio v c /vj as calculated from Eq. (p2[). 
We see that even if the mass transport rate of the jet 
is equal to the mass accretion rate a termination shock 
should form about 50i? o ~ 5 • 10 15 m away from the core 
of the source. This distance is reached by the jet material 
traveling at Uj = 0.61 c in less than two years. 

GRS 1915+105 was discovered as a bright X-ray source 
on the 15th of August 1992 (Castro-Tirado et al. 1992). 
Given the availability of X-ray monitoring satellites be- 
fore 1992 it is not likely that this source was very active 
before this date. The subsequent radio monitoring with 
the Green Bank Interferometer (e.g. Foster et al. 1996) 
shows that after its discovery GRS 1915+105 produced 
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radio outbursts every few months. In the frame of the in- 
ternal shock model described here this implies that jet 
production must have been reasonably steady since 1992. 
Assuming that the bulk velocity of the jet material did 
not vary strongly, the end of the jet must have reached a 
distance of roughly 3 • 10 16 m from the core by the end 
of 1997. This is the time of the radio observation of the 
large scale surroundings of GRS 1915+105 by Rodriguez 
& Mirabel (1998) who did not find any evidence for a ter- 
mination shock of the jet or radio lobes (but see Levinson 
& Blandford 1996a). 

The estimation of the position of the termination shock 
depends crucially on the overdensity of the jet material 
with respect to the ISM. It is possible that the gas den- 
sity in the vicinity of GRS 1915+105 is lower than as- 
sumed here. However, given its location in the galactic 
plane this is rather unlikely. A further possibility is that 
the jets of microquasars are not conical for their entire 
length. The jets of extragalactic FRII objects are believed 
to pass through a very oblique reconfincmcnt shock which 
brings them into pressure equilibrium with their environ- 
ment (e.g. Falle 1991). These shocks are not very efficient 
in accelerating relativistic particles and so are often unde- 
tectable. This scenario is also confirmed for FRII sources 
by numerical simulations of their jets (e.g. Komissarov & 
Falle 1998). The same process may recollimate the jets of 
microquasars as well. In this case they may stay overdense 
with respect to the ISM much longer and this would en- 
able them to travel out to much larger distances before 
terminating in a strong shock. In this respect it is inter- 
esting to note that Rodriguez & Mirabel (1998) found a 
compact non-thermal emission region located 16.3' away 
from GRS 1915+105. The feature is elongated and its ma- 
jor axis is aligned with one of the jets. If this feature is 
caused by the jet pointing in its direction then it must 
have been ejected by the core roughly 280 years ago. This 
may be the time scale on which GRS 1915+105 becomes 
active and produces jets. 

8. Observational tests of the model 

8.1. Emission lines 

In our model the jet flow initially consists of non- 
relativistic hot plasma. Due to adiabatic expansion losses 
of thermal energy during the propagation of this mate- 
rial along the conical jet with given opening angle, the 
temperature of the plasma decreases. Internal shocks as 
envisioned above lead to local heating and acceleration of 
relativistic particles but do not change this general pic- 
ture. In reality the situation is very similar to that in the 
well-known source SS433. In SS433 we observe bright op- 
tical recombination Balmcr, Paschen and Brackett lines 
of hydrogen which are blue-shifted in the approaching jet 
and red-shifted in the receding jet. The velocity of the 
jet flow in SS433 is equal to 0.26 c. Due to the preces- 
sion of the jet the observed red and blue shifts are strong 



functions of time (Margon 1984). In our case the inclina- 
tion angle of the jets of GRS 1915+105 is known from the 
observations by Mirabel & Rodriguez (1994). This angle, 
9 = 68°, and the bulk velocity of 0.6 c of the jet material 
found in our fiducial model permits us to estimate the red 
and blue line shifts of the emitting jet material: 

A ,, , „ „n f 0.97 ; approaching 
-= 7j (l±/3,c OS 0) = | 153 ; re P c P ed . ng (25) 

Note that any line emission coming from the jet approach- 
ing the observer is hardly shifted in wavelength at all. As- 
suming that the bulk velocity of the jet material in the 
jets of GRO J1655-40 is also close to 0.6 c, we find for 
this source that the emission lines are redshifted for the 
approaching jet (A/A' ~ 1.18) as well as for the receding 
jet (A/ A' ~ 1.32). This is caused by the large viewing an- 
gle, 9 — 85°, of the jets in this object (Hjellming & Rupen 
1995). In both cases the very large inclination angles re- 
sult in a strong predicted asymmetry in the line shifts for 
the two jets. Measuring these shifts will permit us to esti- 
mate both the velocity of the jet bulk flow and the viewing 
angle 9. Furthermore, any jet precession as in the case of 
SS433 could be detected. 

Measuring the predicted line shifts is complicated by 
the low density of the material in the jet flow of SS433 and 
GRS 1915+105 which prevents the production of bright 
recombination lines. However, we know that in the case 
of SS433 there is a strong thermal instability in the flow 
which leads to the formation of small, dense cloudlets 
(Panfcrov & Fabrika 1997). This increases the recombi- 
nation rate and effective emission measure of the plasma 
in the flow. If there is a similar instability in the jets 
of GRS 1915+105 and GRO J1655-40 we have a good 
chance to observe recombination lines from both of these 
sources. Another problem is the strong obscuration of 
GRS 1915+105 by interstellar dust. Therefore, it is only 
possible to look for recombination lines of hydrogen in the 
K-band. In the case of GRO J1655-40 obscuration is low 
and there is a chance to detect Lyman and Balmer lines. 
Unfortunately, the mechanical power of the jet in GRO 
J1655-40 is smaller than in GRS 1915+105 or SS433. This 
will lead to a smaller density and emission measure of the 
jet material and therefore also a smaller intensity of the 
lines. It is important to bear in mind that in SS433 the 
emission lines are extraordinarily bright and but modern 
observational techniques permit us to look for blue and 
red-shifted lines which are weaker by many orders of mag- 
nitude. 

In SS433 ASCA discovered red and blue-shifted X-ray 
K- lines of iron with a rest energy of roughly 6.7 keV and 
similar lines of hydrogen- and helium- like sulphur and ar- 
gon (Kotani et al. 1997). In our case the cooling jet flow 
with an initially very high temperature must lead to the 
emission in similar lines of recombining high-Z ions. Again, 
the mechanical energy of the flows in GRS 1915+105 and 
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GRO J1655-40 is smaller than in SS433 and, therefore, 
the lines should be weaker in these objects. However, 
the new X-ray spacecraft, XMM, CHANDRA, ASTRO- 
E, Constellation-X and XEUS, may be able to detect such 
emission in red and blue-shifted X-ray lines. Note in this 
respect the detection of shifted iron lines in GRO J1655- 
40 reported by Balucihska-Church & Church (2000) with 
RXTE which the authors attribute to the accretion disk 
but may very well originate in the continuous jets of this 
source. 

All predictions for the production of line emission in 
the jets are based on the assumption that the jet flow con- 
sists of matter with a high but non-relativistic tempera- 
ture moving as a whole with relativistic bulk velocities. 
There are two other obvious possibilities: (i) The jet mat- 
ter consists of a pair plasma and (ii) the jets consist only 
of ultra- relativistic plasma with no cold electrons present. 
In case (i) we have to consider the possibility of a bub- 
ble around an X-ray source filled with a huge amount of 
positrons. If these positrons become non-relativistic due to 
adiabatic or other energy losses inside the jets and they 
cool down to sufficiently low temperatures, we may ob- 
serve a blue and red-shifted recombination line of positro- 
nium in the optical and UV wavebands. This line has a 
wavelength twice that of the Ly-a line of hydrogen. Much 
more important in this case, annihilation lines could be 
observed again red and blue-shifted relative to the rest 
energy of 511 keV. The strong red-shift but weak blue- 
shift of this line predicted by our model leaves a unique 
signature which will be observable with INTEGRAL. A 
luminosity only a few times smaller than the mechanical 
power of the jets will be emitted in the electron-positron 
annihilation line in this case. This large luminosity should 
make the annihilation lines observable despite the unfa- 
vorable angle of the jets to our line of sight. In the case 
of only ultra-relativistic plasma in the jets, case (ii), no 
recombination or annihilation lines should be observable. 

8.2. Radio continuum 

The internal shock models of GRBs (Rees & Meszaros 
1994) attribute the formation of the shock traveling along 
the jet to the collision of shells of jet material with differ- 
ent bulk velocities. In the non-relativistic limit the velocity 
of the resulting shock is governed by the same momentum 
balance, Eq. (p2|), as the velocity of the termination sur- 
face of the jet. All quantities in that equation with 
subscript 'c' now refer to the slower jet material in 
front of the jet shock while those with subscript 'j' 
denote properties of the faster jet material driving 
the shock. The density ratio r\ is now simply given 
by the densities of the faster jet material driving 
the shock, p2, and that of the slower gas in front of 
the shock, p%. We already pointed out in Sect. 7.1 
that the jet shock is likely to be strong and so both 
Mach numbers in Eq. 22 are significantly greater 



than unity. Therefore v s ~ Vj (l + l/y/rf\ 1 . Since in 
our model all material is assumed to be part of 
the conical jet structure, we find p\ oc p2 oc R~ 2 
and therefore r\ = P1/P2 = const.. This implies that 
within the limitations of the model presented here 
the velocity of the shock is constant as well which 
is confirmed by the observations (e.g. Mirabel & 
Rodriguez 1999 and references therein). 

Once the energy of the shell collision is spent, the shock 
emission fades rapidly. It is therefore possible that we can 
observe the shock reaching the end of the jet only in special 
cases (XTE J1748-288, Hjellming et al. 1999; see above). 
We would then expect that the superluminal component 
should brighten, as well as decelerate rather abruptly. In 
the plasmon model the observed constant superluminal 
motion is taken to indicate a large mass and consequently 
large kinetic energy of the plasmon. If a plasmon is ob- 
served to slow down because of the growing mass of ISM 
it sweeps up, then this deceleration should be rather grad- 
ual unless the plasmon encounters a local overdensity in 
the ISM. The observed deceleration of the superluminal 
component in XTE J1748-288 occurred rather rapidly at 
a distance of about 1" from the core after a phase of ex- 
pansion with practically constant velocity. Furthermore, 
the emission region is still detected in recent observations; 
15 months after the start of the burst (Rupen, private 
communication). During this time it appears to have ad- 
vanced only slowly at a velocity of about 0.01" per month 
or roughly 5000 km s . This slow motion and persis- 
tent radio emission may be interpreted as arising from 
the shock at the end of a continuous jet (see the previous 
section). 

Some interesting predictions can be made from the 
model for future radio observations in the case that these 
can resolve the approaching and receding jet components 
along the jet axis. Because of the way in which the rate of 
acceleration of relativistic particles in the jet by the shock 
varies with time, the peak of the radio emission is not coin- 
cident with the position of the shock. This off-set depends 
on the observing frequency in the sense that the lower 
this frequency the more the emission peak lags behind the 
leading shock. This is illustrated in Fig. |4| where we plot 
the distance of the emission peak on the approaching jet 
side as a function of time for two different frequencies. 
Note also that this effect predicts that we should measure 
slightly different advance velocities of the emission peaks 
at different frequencies. This will not be observed in the 
case of discrete plasmon ejections. 

The steepening of the radio spectrum of the jets in 
microquasars in this model is explained by the superposi- 
tion of the contribution to the total emission from various 
regions within the jet. In resolved radio maps of the jet 
components this should be visible because the model pre- 
dicts the radio spectral index to change along the jet axis. 
This behaviour is shown in Fig. || for the approaching 
jet. Fig. |H also shows the distribution of the flux 
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Fig. 4. Position of the radio emission peak in the fiducial 
model as a function of time after the start of the outburst. 
Solid line: Velocity of the jet shock derived from the 8.4 
GHz observations, dashed line: Model prediction at 1.4 
GHz and long dashed line: Model prediction at 8.4 GHz. 
Only the approaching jet component is plotted. 



along the jet axis. The relatively uniform distri- 
bution is caused by the decrease of the magnetic 
field strength further out along the jet counteract- 
ing the injection of newly accelerated particles by 
the shock. The jet region over which the spectrum 
steepens is small and the emission originating in 
this region also weakens considerably in the di- 
rection away from the jet shock. This may make 
a detection of the spectral steepening along the 
jet difficult. However, the apparent shortening of 
the emission region along the jet at higher observ- 
ing frequencies may be detectable. The decrease in 
the strength of the magnetic field combined with the high 
energy cut-off of the energy spectrum of the relativistic 
particles leads to an overall steepening with time of the 
radio spectrum along the jet axis. This is also shown in 
Fig. H and should be observable if the jet components can 
be resolved at more than one frequency. 

Note also that the length of the region along the jet 
axis which is emitting radiation at a given frequency in- 
creases with time. Although the fraction of the distance of 
the shock from the source centre subtended by the emit- 
ting region shrinks for later times (see Fig. ||), the absolute 
extent of this region will grow. This may also be detectable 
in future radio observations of sufficient surface brightness 
sensitivity. 

Another prediction of the model is that the lightcurves 
of radio outbursts in microquasars at a given observing fre- 
quency should have a fairly constant slope for a few tens 
of days. After that they steepen rapidly once the criti- 
cal frequency of the most energetic relativistic particles 
moves below the observing frequency. This steepening oc- 
curs earlier at higher frequencies. At the same time that 
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Fig. 5. Predicted variation of the radio spectral index be- 
tween 1.4 GHz and 8.4 GHz and flux along the jet axis. 
The length scale on the x-axis is normalised to the dis- 
tance of the jet shock from the source centre at the re- 
spective observing time. Solid line: At the time of the first 
VLA observations, i.e. 4.8 days after the start of the out- 
burst. Short-dashed line: At a time of 40 days after the 
start of the outburst. Dot-dashed line: Flux per relative 
distance in arbitrary units at 8.4 GHz at the time of the 
first VLA observations, corresponding to solid line. Long- 
dashed line: Flux per relative distance 40 days after the 
start of the outburst, corresponding to short-dashed line. 



the steepening of the lightcurves occurs, the flux ratio of 
the approaching and receding jet components should de- 
crease. This effect is seen in Fig. [jj The jet components 
of the second outburst of April 21 observed during the 
second observing campaign with the VLA (Rodriguez & 
Mirabel 1999) show a much steeper lightcurve than those 
of the first outburst and, at the same time, a smaller flux 
ratio. Since these components were observed later in their 
evolution than those of the first outburst, this is in agree- 
ment with the predictions of the model. 



9. Conclusions 

We developed a continuous jet model for the radio out- 
bursts of galactic microquasars. The model naturally ex- 
plains the observed rather flat decaying lightcurves of 
these bursts as the signature of synchrotron radiation 
of relativistic particles accelerated by internal shocks in 
the conical jets. The comparatively long duration of the 
bursts implies that this model is a 'time-resolved' version 
of the internal shock model proposed for GRB (Rees & 
Meszaros 1994), though the synchrotron emission is pro- 
duced at much lower frequencies. The gradual steepening 
of the radio spectrum is explained by a superposition of 
the radiation of different populations of relativistic parti- 
cles with different ages. This spectrum of ages results from 
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the shock traveling along the jet with older populations of 
accelerated particles left behind. 

We find that only a roughly constant rate of accel- 
eration of relativistic particles followed by an exponen- 
tial decay can explain the observed light curves for the 
strong outburst of GRS 1915+105 in 1994. We interpret 
this behaviour as the signature of two colliding shells of 
jet material, as in the internal shock model for GRB. 
A consequence of this is the continued steepening of the 
lightcurves of a given outburst coupled with a decreasing 
flux ratio of the emission observed from the approaching 
to that from the receding jet side. 

The energy requirements of the continuous jet model 
for producing radio outbursts are similar to those of the 
plasmon model. However, much of the energy underlying 
the outbursts may be stored in the continuous jet while the 
passage of the internal shock only 'lights up' the jet. This 
implies that the rate at which the energy of the outburst 
is supplied by the central engine to the jet is much lower 
than in the plasmon model. 

The occurrence of mini-bursts in microquasars with 
flat spectra up to infrared frequencies (Mirabel et al. 1998) 
and the observation of K-band emission in the jet a consid- 
erable distance away from the core (Sams et al. 1996) sug- 
gest different modes of jet production: (i) A stable 'mini- 
burst' mode with relative little variation in the bulk jet 
speed and therefore also only weak internal shocks, (ii) 
A more variable outburst mode with strong variations in 
the jet speed and strong internal shocks (see also Fender 
1999). The weaker flavour internal shocks seem to produce 
flatter relativistic particle spectra extending to higher en- 
ergy compared to the strong shocks. However, the total 
number of accelerated particles must be much larger in 
the strongly variable phase. 

We show that the properties of the continuous jets of 
microquasars should lead to strong shocks at their ends 
where they are in contact with the surrounding ISM. This 
is consistent with the recent observations of the deceler- 
ating radio emission region of XTE J1748-288 (Hjcllming 
et al. 1999) and its persistence for 15 months after the 
start of the original outburst (Rupen, private communica- 
tion). The shocked jet material may subsequently inflate 
a low density cavity around the jets similar to the ra- 
dio lobes in extragalactic jet sources of type FRII. This 
is obs erved in SS433 (Dubncr et al. 1998) and may be 



mon model it would be necessary to spatially resolve the 
super luminal emission regions during outbursts, preferen- 
tially at more than one radio frequency. Additional sup- 
port for the scenario of continuous jets may come from 
further high resolution observations of the cores of mi- 
croquasars during quiescence. These should show at least 
some spatial extension of the radio emission along the jet 
axis as observed in Cygnus X-l (Fender et al. 2000). These 
observations can potentially provide us with valuable in- 
formation on the properties of the jets which otherwise 
we can only study during strong outbursts when strong 
shocks pass through them. 
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